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Qh! Abstract 

. The Garman-Klass unbiased estimator of the variance per unit time 

OC ' of a zero-drift Brownian Motion B, based on financial data that reports 

for time windows of equal length the open (OPEN), minimum [MIN), 
maximum [MAX) and close [CLOSE) values, is quadratic in the statistic 
Oh ■ 5*1 = {CLOSE - OPEN, OPEN ~ MIN, MAX - OPEN) . This estimator, 

•^ ' with efficiency 7.4 with respect to the classical estimator {CLOSE— OPEN)^, 

is widely believed to be of minimal variance. The current report disproves 
2 ' this belief by exhibiting an unbiased estimator with slightly but strictly 

C/3 . higher efficiency 7.7322. The essence of the improvement lies in the proposal 

that the data should be compressed to the statistic S2 defined on W{t) = 
^ . B{0) + [B{t)-B{0)]sigii[{B{l)-B{0)] as S'l was defined on the Brownian path 

K^ ' B{t). The best 5'2-based quadratic unbiased estimator is presented explicitly. 

CN , The Cramer-Rao upper bound for the efficiency of unbiased estimators is 

8.471. It corresponds to the large-sample efficiency of Maximum Likelihood 
estimators. This bound cannot be attained because the distribution is not 
of exponential type. 

f— ^ ' Regression- fitted quadratic functions of S2 (with mean 1) markedly out- 

00 ! perform those of S'l when applied to random walks with heavy-tail-distributed 

^D • increments. Performanceisempirically studied in terms of the tail parameter. 
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1 Introduction 

As stressed repeatedly (see Magdon-Ismail & Atiya (2001)), volatility estimators of financial 
data ought to have as small a variance as possible, because volatilities change over time, 
so past data have decaying importance. The celebrated Garman-Klass (1980) variance 
estimator, introduced almost three decades ago, achieves better accuracy in estimating o"^ 
than the classical, natural estimator average {CLOSE — OPEN)"^ does in seven times the 
observation period. This unbiased variance estimator is the minimum-variance unbiased 
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quadratic function of the spreads c = CLOSE - OPEN, h = MAX - OPEN, I = MIN - 
OPEN (for dose, high, low). These data Si = {c,h,l) can be compressed without loss of 
sufficiency. 

A coarser (but incomplete) sufficient statistic. Consider the triple S2 = (C, H, L) 
where C = \c\ , {H,L) = {h,l) if c > 0, while {H,L) = -{l,h) if c < 0. Without loss of 
relevant information about the variance, the Brownian Motion trajectory {B{t) ; t € (0, 1)} 
may be replaced by the flipped path {W{t) ; t G (0, 1)}, defined as W{t) = B(0) + [B{t) - 
B{0)]sign{B{l) - B{0)). That is, the three interval lengths {-L,C,H - C), in fact the 
further compression (C, min(— L,/:/^ — C),max(— L,i:/^ — C)), determined by {c,h,l), carryall 
relevant information contained in (c, h, I) about o"^, but do not determine (c, h, I). Although 
intuitively clear after some thought, sufficiency of (C,mm{—L,H — C),inax{—L,H — C)) 
can be formally inferred from Siegmund's (1985) representation displayed as (14) in the 
sequel. The Rao-Blackwell theorem (Blackwell (1947), Rao (1946)) claims that under these 
conditions, for every S'l-based unbiased estimator of some arbitrary parameter there is 
an S'2-based unbiased estimator with smaller variance - strictly smaller unless the two 
coincide. As will be seen, the Garman-Klass estimator is a function of 5*2, so the Rao- 
Blackwell improvement leaves it invariant. However, the Garman-Klass estimator, best 
among the quadratic function of Si, is not best possible as a function of 5*2. Had 5*2 been 
a complete minimal sufficient statistic, Garman-Klass and the proposed estimator would 
have equally been the UMVUE (uniformly minimum variance unbiased estimator) of the 
parameter. However, (7^ and 2[{H — C)"^ + L^] are different unbiased estimators of o"^. 
Hence, 5*2 (whether minimal sufficient or not) is not complete. Loose some, win some: 
we will only conjecture rather than claim optimality of the proposed S'2-based quadratic 
unbiased estimator of o"^; on the other hand, the exchangeability property under which 
(— L, C,H — C) and {H — C, C, —L) are identically distributed, justifies searching for the 
best quadratic function of {—L, C,H — C) among those that are linear combinations of four 
rather than six quadratic terms. 
Four basic quadratic unbiased variance estimators. 



ai = 2[{H-CY + L'], ai = C\ di = 2{B - C - L)C , ai = -^-—-^ (1) 

21og(2)-| 

The rationale for the somewhat bizarre coefficients is that each of these four terms is an 
unbiased estimator of o"^ , with respective variances 

Mw{a\) = 0.797943^^ , Var(a|) = 2o-^ Nw{al) = 0.504753ct^ , Var(a|) = 1.004876^^ (2) 



The proposed variance estimator vis a vis Garman Klass. The proposed estimator 
(T^ = X^i <^i^i assigns to these four terms respective weights 

ai = 0.273520 , 02 = 0.160358 , 03 = 0.365212 , 04 = 0.200910 (3) 

and achieves variance Var((T^) = 0.258658cr'^. The Garman-Klass estimator [3] 

6-2,^ = 0.511(/i - If - 0.019(c(/i + /) - 2hl) - 0.3830^ (4) 

happens to pool these four basic estimators too, so the Rao-Blackwell theorem does not 
rule out the possibihty that it coincides with a^. However, as argued earher, the two do 
not agree, and a^^ = J2i Pi^^f pays a price for being quadratic in {c,h,l). Its coefficients 
are given by 

0.511 
Pi = — — = 0.2555 

(32 = 0.511 - 0.383 - 0.019 = 0.1090 

0.019 
p3 = 0.511 — = 0.5015 

/?4 = 2(0.511 -0.019)(21og(2) --) = 0.1340 (5) 

that achieve Var((T^^) = 0.27o"^. 

Maximum Likelihood variance estimators and Fisher information. In principle, 
giving up on the requirement of unbiasedness, the computer-intensive maximum likelihood 
estimator (MLE) of cr^ by Magdon- Ismail & Atiya (2001) could have been a competitor, since 
MLE's are functions of any sufficient statistic. However, this estimator is based on (h, I) 
rather than on (c,h,l). Magdon- Ismail & Atiya report that their estimator has variance 
slightly higher than Garman-Klass'. 

The joint generating function of (c, h, I) is presented by Garman & Klass as an infinite 
series, from which these authors derived all pertinent second and fourth degree moments. 

Ball &: Torous (1984) developed an infinite-series formula for the joint density of (c, h, I) 
and used it to construct numerically the MLE of o"^. They report estimated efficiency of the 
MLE for a selection of sample sizes, basing each value on a simulation sample size of 1000 
runs, a great achievement in 1984, but insufficient for delicate comparisons. An attempt 
at numerical evaluation of the Fisher information, based on the Ball &: Torous expression 
for the joint density, disclosed that their formula seems to have a missprint. This joint 
density was re-derived based on the formula by Siegmund quoted earlier, exhibited as (14) 



in the sequel. The inverse of the Fisher information is the Cramer-Rao lower bound for 
the variance per time-window of any unbiased estimator of o"^, for any sample size. It is 
also the asymptotic variance of the (not necessarily unbiased) MLE of a^. Its value turns 
out to be 0.2361. This is the benchmark with which Garman-Klass' 0.27 and the proposed 
estimate's 0.258658 variances should be compared. 

The Cramer Rao bound 0.2361 is not attained by unbiased variance estimators: 
disproving exponentiality of a family of distributions. Under proper regularity 
assumptions (see Joshi (1976) ), the Cramer-Rao bound is attained if and only if there is 
a linear relationship between the estimator and the score function (derivative with respect 
to the parameter of the logarithm of the density). However, for this to happen, there must 
exist a linear relationship between the score functions evaluated at different values of the 
parameter. It was ascertained numerically that this is not the case. In other words, the 
model is not of exponential type. We don't know whether the sufficient statistic 5*2, shown 
above not to be complete, is minimal sufficient. As a result of all of these considerations, 
the proposed estimator may not be of minimal variance. 

Since both the proposed and Garman-Klass' estimators are averages over time-windows, 
their variances per time-window are independent of sample size. It is conceivable, and Ball 
& Torous have provided evidence in this direction, that the MLE has variance per time- 
window that decreases as the sample size increases, so for small sample sizes the proposed 
estimator has in practice no competitor. 

Moreover, since the BM model doesn't really hold in practice, a broader contribution of 
this paper is the introduction of more efficient quadratic statistics on which to base practical 
estimators. Simulation results for random walks with t-distributed increments are reported 
in Section 3. 

2 Derivation 

Following the steps of Carman & Klass, all second and fourth order moments of (C, L, H) 
will be identified. Some of these will be quoted from Carman &: Klass, some will be derived 
once the joint densities of (C, H) and (C, L) are explicitly presented, and some will require 
some additional argument. Although it would perhaps be more natural to work only with 
the exchangeable variables A = H — C and 6 = —L, work will be performed on the variables 
H and L as well, in order to link more easily with Carman &: Klass' triple (c, h, I). 



2.1 The joint densities of C and each of H and L: four unbiased estimators 

Assume throughout the computations that the drift is and the variance per unit time is 
1. Thus, E[C^] = E[c^] = 1. 

By a common reflection argument, BM reaching at least as high as a; > and ending up 
at y = X — {x — y) G (0,x) is tantamount to ending up at x + (x — y). Or, 
P{H > x,C G [y,y + dy]) = P{C G [2x - y,2x - y + dy\) = 2(t){2x - y)dy, where 
(f){-) = —7= exp{— 2(-)^} is the standard normal density function (see Siegmund or expression 
(14) in the Appendix for a generalization to (C,H,L)). 

Similarly, P{L < z,C £ [y,y + dy]) = P{C G [2z-y,2z-y + dy\) = 2(j){2z - y). Hence, 
the joint density of H and C is 

fH,c{x, y) = 4(2x - y)(t){2x - y) , < y < x (6) 

and that of L and C is 

fL,c{z, y) = 4(y - 2z)(t>{y - 2z) , z < < y (7) 

These joint densities, essentially re-phrasings of a well known formula for the joint 
density of (ft,, /i — c) (see Yor (1997)), lead to the flrst four of the following five second 
moments. The fifth is taken from Garman & Klass. Details are omitted. E[C'^]=1 by 
assumption. 

E[H^] = I , EXL^X = \ , E\Cli\ = \ , E[CL] = -\ , E[HL] = 1 - 21og(2) (8) 

As a corollary, 

Lemma 1 The variance estimators di , i = 1,2,3,4 (see {!)) are unbiased. 

Seshadri's (1988) theorem that 2h{h — c) is exponentially distributed with mean 1, and 
is independent of c, implies that 2H{H — C) is exponentially distributed with mean 1, and 
is independent of C. This is so, simply because the conditional distribution of (/i, c) given 
that c > is the (unconditional) distribution of [H^ C). 

Of course, the same applies to 2l{l — c) and 2L{L — C). However, 2H{H — C) and 
2L{L - C) are dependent (identities (10) yield correlation 1 + |C(3) - 81og(2) = -0.3380 
between the two), and dependent given C. 

Otherwise, it would have been very easy to sample (C, H, L) triples. As things stand, it 
is easy to sample pairs (c, h) (and (c, /)) or (C, H) (and (C, L)), by independently sampling 



c and h{h — c). A practical approximate method to sample {C,H,L) triples is to sample 
{C',H') correctly, then make the wrong choice L' = C — H' , not on [0, 1] but on each of 
the N sub-intervals [^^,j^]- The construction is correct except if H and L are attained 
in the same sub-interval, the probability of which decreases fast as N increases. Instead of 
letting L' = C — H' , other copulas may be used, to better approximate features of the joint 
distribution of (C",i7',L'). 

2.2 The MLE's of a^ based on (C, H) and on (C, L) are unbiased 

It may be of interest to notice that (6) (resp. (7)), reinterpreted as fH,c{xiy'-,o') = 
4^fi^</>(^^), identifies the MLE of cj^ based on {C,H) (resp. (C,L)) as the average 
over the sample of \{2H - Cf = ^C^ + \[A{H - Cf] + \[^C{H - C)] and \{2L - Cf = 
oC^ -)- 3[4L^] -|- g[— 4C-L]. The average of the two, the simple average of the first three 
unbiased estimators in (1), achieves variance 0.3694, above Garman-Klass'. 

2.3 The fourth moments of (C, H, L) 

The following fourth moments are derived from the joint densities of {H,C) and (L,C). 
ElC^] = 3 is Gaussian kurtosis. 

27 3 31 1 

E[C'H] = - , E[C'L] = -- , E[C'H'] = - , E[C'l'] = - (9) 

The following fourth moment information is taken from Carman & Klass. (^ is Riemann's 
zeta function, with C(3) = EfcLi pr ~ 1.2020569. 

E[H^L'^] = ^[/i2?2] =3_41og(2) 

E[C^HL] = E[c^ hi] =2 -2log{2)-lc{3) 

8 

E[H^L] + E[HL^] = E[hl{h^ + f)] =6-6 log(2) - ^C(3) 
E[CH^L] + E[CHL^] = E[chl{h + 0] = | - 41og(2) - Jc(3) (10) 

There is one more (C, H, L)-based fourth moment needed, whose value does not follow 
from Garman &: Klass'. 

Lemma 2 E[CHL^] = C(3)/16 - 21og(2) + g ~ 0.1575842. 



A proof of Lemma 2 can be found in the Appendix. Large sample empirical estimation 
of E[CHL?] gave 0.15762, yielding \ai{a1) very close to 1. Had E[CHL?] been equal to 
log(2)(3 — 41og(2)) ~ 0.15763 (initial conjecture), Var((T4) would have been exactly 1. 

From all the fourth moments above, 



E[C^] 


= 


3 


E[5^] 


= 


^[.^]^^ 


E[C5^] 


= 


-^[c.3].A 


E[C^5^] 


= 


E[C'L^] = \ 


E[C^5] 


= 


-E[C'L] = I 


E[C^M] 


= 


E[C'L] - E[C^HL] = 21og(2) + Jc(3) - 

o 


E[CM^] 


= 


E[CHL^] - E[C^L^] = E[CHL^] - I 

8 




= 


43 
C(3)/16-21og(2) + - 


El^^d"^] 


= 


E[H^L^] + E[C'^L^] - 2E[CHL^] 






3 C(3) 






16 8 


2E[^^5] 


= 


E[/\^5] + E[/^5^] = -{E[H^L\ + E[HL^ 




+ 


E[C^L] + E[CL^] - 3E[C^HL] + 3E[Cl 




= 


61og(2)-A^(3)-^ 



19 



(11) 



2.4 The covariance matrix of the four basic estimators 



Let S stand for the covariance matrix of the four basic estimators. Their variances are on 
the diagonal, their covariances off the diagonal. 

Applying the formulas of the previous sub-section, the variances of the basic estimators 
af (see (1)) are 

S(l,l) =Var(af) = 8(^[(5^] + ^[A2<52]) _ i = 2 - C(3) = 0.797943 

S(2,2) = Var(a|) = 3-1 = 2 

S(3,3)=Var(ai) = 8iE[CH^] + E[C^A5]) - 1 = 8[logi4) + '^-((3) - ^] - 1 



= 0.504753 

pr a2 r2i .3 C(3) 

E(4,4)=Var(^|) = , J,, _ L - 1 = ,"_%,, - 1 = 1-004876 (12) 



The covariances of the basic estimators are 

S(l,2) = Cov(<T?,<Ti) = AE[CH^]-l = -]^ 

S(l,3) = Cov(a?,<Ti) = 8^[C53]+8E[CA52]-l = ^ii-^-161og(2) 

= 0.010674 

AE[^5^] ^_121og(2)-f -|C(3) ^ 



2 ^2\ 



S(l,4) = Cov(af,(Ti 



log(4) - I log(4) - f 

.580786 



S(2,3) = Cov(ai,ai) = AE[CH]-l = ]^ 



log(4) - I log(4) 4 



S(2,4) = Cov(ai,a|) = -^H^-^ - 1 = ^1^^ = -.537074 



4.E[C^^S\ _ _ ^ + f - 81og(2) 
log(4)-|" " log(4) 4 



S(3,4) = Cov(ai,^l) = ^^^=,^-1= ' , ', r' -^ 



= -.043711 (13) 

2.5 Derivation of the proposed estimator 

Letting a (see (3)) stand for the weights assigned to the basic estimators, the weighted sum 
has variance a^Sa and mean a^l. Using a Lagrange multiplier to constrain the mean to 
be 1, minimal variance is achieved at a = iTy;-ii ; yielding the weights displayed in (3). 
The variance of the proposed estimator is ^t^-i-^ = 0.258658, with corresponding efficiency 
21'^S^il = 7.73221. 

3 Heavy tailed random walks - simulation results 

As is commonly observed in financial data, the logarithmic increments of returns have 
power-law tails, at least in the visible range, with tail parameter around 3. This means 
finite variance but infinite variance of the usual empirical variance estimators. Suppose 
that the basic process on which (Open, Close, Min, Max) data is reported per time window 
is a random walk with t-distributed increments. A simulation analysis will now be reported, 
in which the number of increments of the random walk per time window is 10, 30 and 50, and 
the degrees of freedom (df) range from 1.5 to 5 with step size 0.5. Minimum sum-of-squares 
quadratic functions with mean 1 of the 5*1 and 5*2 statistics were fitted by Regression, with 
sample size 10^: the regression coefficients were identically calibrated so that the predictor 
of unity has mean 1 in each such sample. Each such Regression was repeated 100 times, 
and the averages of the corresponding regression coefficients and overall "variances" were 



recorded. Of course, second moments are finite only for df > 2 and fourth moments are 
finite only for df > 4, but the empirical study seems instructive. A sample of size 10^ from 
the sum of iV = 50 t|(i/-=3|-distributed random variables typically displays lighter tails than 
df = 3 would entail. Table 1 reports the empirical minimum variance of the quadratic 
functions, and Table 2 reports the coefficients of the building blocks of expression (1) that 
yield the minimum-variance quadratic function for each case. These building blocks have 
expectation 1 for Brownian Motion but not for random walk, so their coefficients need not 
add up to unity. Table 1 displays performances similar to those derived for Brownian Motion 
for moderate df, fast deteriorating when df decreases, in which case 5*2 data progressively 
outperforms 5i data. S2 data yields lower variances than Si data throughout the range, as 
well as for uniform and double exponentially distributed increments, although the difference 
in variance in these light-tail cases is as small as for BM. 



Table 1. Minimum variance of mean-1 quadratic functions of Si and 5*2 data 



df N 


10, ^2 


10, Si 


30, ^2 


30, Si 


50, ^2 


50, Si 


1.5 


16.2403 


51.0366 


8.3438 


32.4697 


6.5322 


28.3950 


2.0 


4.8444 


6.6039 


2.6532 


3.8327 


2.1972 


3.2252 


2.5 


2.5864 


2.8365 


1.4297 


1.5529 


1.1718 


1.2627 


3.0 


1.7359 


1.8038 


0.9527 


0.9782 


0.7630 


0.7788 


3.5 


1.2334 


1.2746 


0.6809 


0.6991 


0.5467 


0.5624 


4.0 


0.9469 


0.9776 


0.5409 


0.5585 


0.4532 


0.4686 


4.5 


0.7864 


0.8124 


0.4792 


0.4957 


0.4094 


0.4239 


5.0 


0.7071 


0.7296 


0.4473 


0.4629 


0.3896 


0.4037 


CX) 


0.4679 


0.4826 


0.3630 


0.3765 


0.3369 


0.3496 


00, A^ = 00 










0.2587 


0.27 



It is of interest to observe how does 6*2 outperform 5*1 data for low df. Table 2 shows 
that the role of C is downplayed or even dampened in favor of those oi H — C and —L, 
gradually incorporating C into the Regression as df increases. The rationale for this is that 
the tail parameter of sums of i.i.d. data is the same as that of the summands, whereas the 
tail parameter of extrema is the sum of those of the summands. This makes C theoretically 
as heavy tailed as each increment, but makes H — C and —L have lighter tails than the 
increments. In contrast, the [h, c, I] data of statistic Si is less able to split variables into 
light tail and heavy tail components. Although h — \c\ — I = H — C — L, the insistence 
on resorting to quadratic functions leaves it out of the Si game. Still, both statistics seem 
to work fairly well even under low df. In contrast to the variances 2.1972 or 3.2252 for 
df = 2, 0.7630 or 0.7788 for df = 3 and 0.4532 or 0.4686 for df = 4 (see A = 50 in Table 



1), the calibrated C^ has respective empirical variance above 5000, 16 and 2.5, converging 



reasonably fast (2 + 



6 

(4f-4)Af. 



to 2 thereafter. 



Table 2. Coefficients of the minimum variance mean-1 quadratic function of 5*2 
data for A^ = 50 increments per time window 



df 


2{{H-Cf + L^) 


C2 


2{H-C-L)C 


-(U-C}L 

21og(2)-5/4 


1.5 


0.0209 


-0.0000 


0.0010 


0.1724 


2.0 


0.1358 


-0.0004 


0.0352 


0.1561 


2.5 


0.1745 


-0.0034 


0.1573 


0.1215 


3.0 


0.1827 


0.0140 


0.2461 


0.1149 


3.5 


0.2006 


0.0666 


0.2460 


0.1228 


4.0 


0.2185 


0.1081 


0.2442 


0.1317 


4.5 


0.2335 


0.1271 


0.2620 


0.1399 


5.0 


0.2480 


0.1395 


0.2781 


0.1473 


oo 


0.3974 


0.2321 


0.4390 


0.2245 


oo, Af = oo 


0.2736 


0.1604 


0.3652 


0.2009 



4 Appendix - proof of Lemma 2 

For the sake of conciseness, the tedious integration to be presented will be restricted to the 
identification of E[CHL'^], although, in principle, more general joint moments and moment 
generating function of (C, H, L) could have been identified. 

Consider the infinitesimal event {BM{1) € (CC + ^O , BM{s) G (a, 6) , Vs G [0,1]}, 
where a < min(^, 0) < < max(^, 0) < b. By Siegmund's Corollary 3.43, its probability 
Q{S,,a,b)d^ is as follows 



Q{ta,b)= J2 {HC-2j(b-a))-cP{C-2a-2j{b-a))} 



(14) 



j=-oo 



The joint density fc,h,iiS,, ci, b) is (minus) the mixed second derivative of Q with respect 
to a and 6, on {^ G (a, 6) , a < ,6 > 0}. The joint density fc,H,L is simply 2fc^h,i, 
restricted to {^ G (0, 5) , a < ,b > 0}. The two terms in the j = and second term in 
the J = 1 summands vanish because they are independent of at least one of a and b. 

To calculate E[CHL'^], the contribution of each summand in (14) will be integrated in 
three univariate steps. The first step will integrate over a G (— oo, 0) the product of a^ and 
the pertinent mixed second derivative. -§^4^{i + Ka + Mb)da is to be interpreted as the 
integration-by-parts element (i<^(^ -|- Ka + Mb), viewed as a function of a. 

'•0 d o d 



I -KTa^^(t){i + Ka + Mb)da 
J-oo ob da 



= i^4['^(? + M^) + (? + Mb)<^{C + Mb)] (for i^ > 0) 

= ^<5(e + M6) (fori^>0) 

= ^HC + Mb)-^{ioTK<0) (15) 

Now expression (15) will be multiplied by S, and integrated over ^ € (0,6). For A > 
(-ftT < 0) it is convenient to integrate ^* {^). These terms appear in (16) and (17). The 
free term in (15) contributes -ttt-j" ^^'^ cancels with the corresponding b^ term in (17). 



f ^T^ f a^ct){i + Kda + Mb)di 

Jo db J-oo 

r(M+i)b 2M% r{M+i}o 

/ y<^{y)dy - —-^ \ ^{y)dy 

J Mb J^ J Mb 



/O OO J-oo 

2M fiM+i)b _ _ 2MH f(M+i)b 

M 

+Mb4>{Mb) - (M - l)b4>{{M + 1)6)] 
M 



— 2 [(M^ft^ + l)$(M6) - ((M^ - 1)^2 + i)$((M + 1)6) 



[(m262 + 1)$* (M6) - ((M^ - 1)6^ + 1)$* ((M + 1)6) (16) 

+Mb(P{Mb) - (M - 1)60((M + 1)6)] + ^6^ (17) 

Finally, expressions (16) and (17), multiplied by 6 and integrated over 6 € (0, oo), via 

yield a rational function of j (with M = 2j and K = —2j or K = —2{j — 1)) whose 
sum contains only terms of the form — X]i°(~l)"'^ = log(2) and J^T ^ = C(3), as in the 
statement of Lemma 2. Further details are omitted. 
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